Photocurrent in nanostructures with asymmetric antidots 
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The steady current induced by electromagnetic field in a 2D system with asymmetric scatterers 
is studied. The scatterers are assumed to be oriented cuts with one diffusive and another specular 
sides. Besides, the existence of isotropic impurity scatterers is assumed. This simple model simulates 
the lattice of half-disk which have been studied numerically recently. The model allows the exact 
solution in the framework of the kinetic equation. The static current response in the second order 
of electric field is obtained. The photogalvanic tensor contains both responses to linear and circular 
polarization of electromagnetic field. The model possesses non-analyticity with regards to the rate 
of impurity scattering. 

PACS numbers: 73.40.-c, 73.50.Bk, 73.50.Pz 
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INTRODUCTION 

Nowadays technology allows to fabricate artificial ar- 
rays of antidot scatterers which form superlattices in 
semiconductor heterostructures with a two dimensional 
electron gas (2DES) (see e.g. Ul l2l ll Lt and Refs. 
therein). The size of antidots can be varied from a few 
microns to a few tens of nanometers at a typical electron 
density n e ~ 10 12 cto~ 2 . Superlattices with circular an- 
tidots (disks), known as the Galton board have been 
realized in experiments 01 II & ■ The mathematical 
theorems of Sinai garanty that the classical dynamics of 
electrons in such structures is chaotic Q. The effects 
of chaotic dynamics and contributions of unstable pe- 
riodic orbits on electron conductivity have been clearly 
seen experimentally. They were also analyzed by theo- 
retical methods and numerical simulations in great detail 
■ The irradiation of such superlattices by a microwave 
field H opens interesting possibilities for microwave con- 
trol of electron current in nanostructures. These systems 
of relatively large period are classical, and the periodic 
potential is the source of electron scattering. 

Unlike arrays of symmetric antidots, studied experi- 
mentally in early works, more sophisticated systems were 
for a long time out of attention. Meanwhile, systems 
without inversion symmetry are capable to rectify the 
electric current. The stationary current in homogeneous 
media affected by light in the absence of any dc-voltage 
called the photogalvanic effect (PGE) was studied since 
the end of 70th 0, ES El 111 E Q • Similar direct cur- 
rent caused by temporal irreversibility due to simultane- 
ous action of two electromagnetic fields with frequencies 
uj and 2lo is known as the coherent photogalvanic effect 
(CPGE)) Q. 

Appearance of directed transport without obvious di- 
rected forces is also known as ratchet effect which has a 
long history. For example, the behavior of a ratchet un- 
der the influence of thermal fluctuations was considered 
in the textbook by Feynman, Leighton and Sands [ljj in 
connection with the problem of reversibility in statisti- 
cal mechanics. Recently the ratchet problem attracted a 



great interest of the scientific community [l7l fisj l . In fact 
the photogalvanic effect corresponds to a ratchet sub- 
jected to weak alternating force with zero mean. Such 
ratchets have been observed in various physical systems 
including Josephson junction arrays 0, |H EL cold 
atoms |22j, macro por ous silicon membranes |23| . mi- 
crofluidic channels [24j and other systems. The growing 
interest to ratchets is strongly stimulated by t heir pos- 
sible applications to biological systems 0, [2^ . In this 
sense the artificial asymmetric nanostructures, as those 
discussed in [H HI Esj and here, can serve as a pro- 
totype for understanding of photocurrent properties in 
biomolecules. The results obtained for ratchets induced 
by microwave fields in nanostructures can be also used for 
understanding of directed transport created by ac-fields 
in molecular electronics [2^|. 

Relatively recently an artificial lattice of asymmetric 
antidots (triangles) participating in transport as scatter- 
ers of electrons was realized experimentally. It exhibited 
the direct current induced by alternating electric field 
I30l . The effect of such type was considered theoretically 
in 31] by means of the classical kinetic equation for the 
system with weak asymmetric periodic lateral potential. 

Another case of superlattices of asymmetric antidots 
(semidisks oriented in one direction or the semidisk Gal- 
ton board) has been proposed and analyzed theoretically 
HE HI El| by simulations of motion of a particle sub- 
jected to alternating force with zero means and colli- 
sions with hard-wall antidots. It has been shown that 
a microwave radiation creates the directed flow of elec- 
trons. The velocity of this flow Vf is proportional to a 
friction coefficient while the direction of current depends 
on the microwave polarization. The directed transport 
induced by a microwave field appears also for an ensem- 
ble of noninteracting particles been at a thermal equilib- 
rium at temperature T. This effect is absent in struc- 
tures with circular antidots due to symmetry conserva- 
tion. Recently, semidisk Galton board realized experi- 
mentally [Hi exhibited the rectification of high-frequency 
electric field. 

In [il |H the properties of photocurrent in asym- 



2 



metric nanostructures have been explained on the basis 
of heuristic arguments and extensive numerical simula- 
tions. However, the analytical theory of the effect still 
needs to be developed. This aim is reached in this paper 
with the help of a kinetic equation approach applied to 
a specific model of asymmetric scatterers (cuts model). 
This approach allows to obtain analytical dependence of 
photocurrent properties on system parameters. 



KINETIC EQUATION APPROACH 

Here we study a simple model of anisotropic 2D arti- 
ficial scatterers, which permits the analytical considera- 
tion in the framework of kinetic equation approximation 
and leads to the PGE. The kinetic equation was used 
in a number of papers devoted to photogalvanic effect. 
Unlike the particle dynamics method, this way suggests 
the developed chaos picture where ergodicity of motion is 
achieved. At the same time it is free from the voluntary 
assumptions about the particle friction used in ^ . 

The system under consideration contains randomly 
distributed oriented scatterers. The scatterers are as- 
sumed to be segments of the length D oriented along the 
y axis; one side (left) of the segment is specular and the 
other side is diffusive (see Fig. P). The concentration 
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approximation the kind of spatial distribution of scatter- 
ers (random or periodic) is of no importance. Besides, 
to limit the possible divergency of the result we include 
isotropic impurity scattering into our model. 

This system can be considered as a simplification of 
the semicircle model studied in j2|| [2]J. In fact, the 
diffusive side of the cut scatters particles like round side 
of the semicircle, if not to pay attention on the difference 
between randomized (in our case) and deterministic (in 
semicircle case) motion. The advantage of " cuts" model 
is its exact solvability. 

The kinetic equation reads as 



the distribution 



(1) 



where f{p,f) is the distribution function, p = 
p(cos if, sin ip) is the electron momentum. The term (Ff) 



represents the action of electric field E(t) = Re(E a 
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of the electromagnetic wave with the complex amplitude 
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(2) 



As it will be seen further, the acceleration in the y di- 
rection can not be limited by the scattering on segments 
only because electrons moving in this direction do not 
relax. This is why we have not restricted our considera- 
tion by the collision with cuts only but have also taken 
into account the impurity scattering. Hence the colli- 
sion integral If = I%f + I c f is assumed to consist of 
impurity-induced 7, and segment-induced I c collision in- 
tegrals. The segment-induced collision integral 



FIG. 1: Considered model system. Cuts have specular left 
sides and diffusive right sides. This produces anisotropy of 
scattering resulting in the photocurrent. 

of scatterers N is supposedly low: ND 2 < 1, In this 
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I c f= / d<p'[W{ip', V )f(p,<p')-W{ V ,<p')f(p,v)}. (3) 
Jo 

is determined by the scattering probability on the cuts 



W(<p',<p) 



cos ip' #(cos (p')S(cp' + ip — 7r) 

r 



■ COS ip' cos ip #(cos (fi)9(— COS ip') 



(4) 



where r = (DNv) -1 is the corresponding characteris- 
tic time, v — p/m; 9(x) is the Heaviside function. The 
impurity scattering is suggested to be isotropic: 



hf 



1 



(/- < / >), 



(5) 



where < ... > means average over angles, Tj is the relax- 
ation time corresponding additional isotropic scattering. 

We shall solve the kinetic equation Eq.QJin the sec- 
ond order on electric field. In the elastic approximation 
the collision operator is degenerate because its action on 
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the isotropic function gives zero; thus this operator has 
an eigenfunction with zero eigenvalue. Physically, this 
means that the elastic scattering does not change any 
distribution function depending only on the energy. But 
if to project the Gilbert space of distribution functions 
on the subspace with zero angular average, the corre- 
sponding projection of the scattering operator becomes 
non-degenerate. So if to consider only the distribution 
functions with zero mean, the elastic scattering opera- 
tor remains non-degenerate and can be treated as full 
relaxation operator. Inelastic scattering controls only 
isotropic part of the distribution function which (if it is 
weaker than elastic one) has no impact on the PGE. 

In the second order on alternating electric field the 
steady current density is described by the phenomeno- 
logical expression 

ji = OtijkE^jE-u,):- (6) 

The symmetry of the considered system allows the fol- 
lowing non-zero components of the photogalvanic tensor 

O-xxxi a xyyi a yxy — a yyx- CO 

From Eq.© it follows that components a xxx and a xyy 
are real. The same is valid for all components of c&ijk in 
the static limit (u> = 0). 

The equation @ is specified as 

jx ^xxx\E x \ "t" OL x y y | Ey | , 

j y = Re(ayxy){E x E; + E*E y ) + Im(a^)[EE*] z . (8) 



The components a XX x, &xyy and Jie(a yxy ) determine the 
response to the linear-polarized light. For linear polar- 
ization along x or y axes the current flows along the x 
direction; the current in y direction appears for tilted 
linear-polarized electric field. In the case of circular po- 
larization the y component of the current is determined 
by \m.(a V xy) (the circular photogalvanic effect) and by 
the sign of the rotation, while x component of the cur- 
rent is the sum of responses to x and y linear polarized 
light and does not depend on the sign of [EE*] Z . No- 
tice that the circular photogalvanic effect vanishes if the 
frequency u> — > 0. 

The formal solution of the Eq. (|TJ in the second order 
in electric field is 



h = {--i)- 1 F(--i)-'Ff, (9) 



where /q = l/(exp (e — /i)/T) _1 is the isotropic equilib- 
rium distribution function, [i is the chemical potential, T 
is the temperature. The kernel of the inverse operator in 
frequency representation, the Green function G w (ip',ip), 
satisfies the equation (— iu> + I)G u (tp, if') — 8 (if — if'). 
This equation can be solved exactly: 



G*W) = -t{ 

0(— cos ip) 
(|cosy|+£) 2 



5{ip — ip') | cosy | 



COS 2 If' 



£ cos ip' 



- COS If 



I cosyl + £ 2 
- cos if 5 (if + if- — 7r) + £ cos if 

I 



(|cosy|+£)(|cosy'|+£) 2 (l -<&(£)) 
cos 2 f' — £6*(— cos if') cos f' 



2(|cosy'|+£) 2 (l -*(0) 



(10) 



The operator G w depends on the frequency via the pa- 
rameter £ = —iiOT + t/Tj, 



$(x) = —TtX + 2 - 



1 



, 3 - 2x 2 



1 — x 2 1 — X 2 



(11) 



The functions denote the integrals expressing via 

elementary functions: 



3>(n, x) 



*' 2 dt 



o (cost + x) n 



arctanh( — ) (12) 



The stationary part of the correction fa can be written 



as 



1 



h = -G°(F U G-"F_ U + F_ u G"F u )f 0] G° = G"=°. 

(13) 

The static current reads 

ji 



, , d 2 p - 



(14) 



The photogalvanic tensor is expressed via partial tensor 
&ijk for electrons with given p : 



2mir'' 



f , a , df Q 

dp t (-e—)a ijk . 



(15) 
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In the degenerate case the partial tensor itself determines The dimensionless tensor aijk can be presented as a func- 

the total tensor a^k- tion of parameters £ = r/r, and £ = £ — ifl (CI = ur). 

aijk(T — 0, cf ) = «o a»jfc|p=pj?) (16) 

where a = e 3 /(47r 2 ^ 2 L> 2 w F ), uf = p F /m,p F is the ANALYTICAL RESULTS 

Fermi momentum, ep = p F /2m is the Fermi energy. At 

finite temperatures T Substituting the stationary part of distribution func- 

/g j. tion (|13() into the expression (|14() and using the Green 

de ( — g— ) Q!jjfc(T = 0, ejr = e). (17) function I|1U|) we find after integration: 



=-^u i^n | i-^ + ^ + A J C % (i,fl- C * 2C ^ 4 4(U) + £ 2 ^-^ 



Bio B (o ( - $ - ^ ^y i ^\^ o + 2 ^^ 4+c2( -^g- c(e+3e3)) .a, o - 

CMC^-2C^) , (2! + ^(,-4^C(-4 + 70) . (2| ^ + ^±±i% (3 , ) + ,, (18) 
= 1 - 35(C) + <( 4j? (fl-l) + (3C 3 ^ - C 4 - + C 2 K 2 i 2 + B(C)(K 4 + 3« - 15C 3 C - 3C 2 + lO^ 2 )) 

1 Is - 

| g(c 2 -ce -2c 2 g 2 + 3q 3 - e 4 + g(0(3ce-3c 2 + ioc 2 e 2 -i5ce 3 + 5^)^(1, e) | c 2 £(0(c 3 + c- 2^)4(2,0 

« - o 3 (e - C) 2 

e 2 (c - g - ce + e 3 + B(cm - k + s^ 2 - ?a) o geg^) (g - 1)) $ (3 , ^ 

«-o 2 c-c 

<W6 o = ^(c - 0(-2« + ^ 2 c(c - 0(-i + e 2 ) + *(-£ + c(2 + (c - 360)) + + 

2C 2 (1 + + 2C 2 ) + C 3 (tt - 2< 2 ) + 2(£(2 + 2< - 3£ 2 - 2< 3 ) + 4C£ 3 (-1 + £ 2 )4(1, 0) + 

C(C - 20$(1, C)("2(-l + C 2 )e - tt(-1 + 2C 2 )(-1 + 0(1 + + 2C(-e 2 + C 2 M + 2e 2 ))$(l,0) ) x 



(19) 



1 



^2(C-0 2 (^ + £-< 2 +£(-3 + 2^(1,0) J -I (j , 

2(c 5 + 2^ - 3c 2 e - 3c 4 e + c(2 - 4e 2 ) + c 3 (-i + 6$ a ) w> + 
2(-(e(-2 + e 2 )) + c 2 e(-2 + 3£ 2 ) + c(2 - be + e 4 ))$(i, o 



(C - 0(-8 + < 3 - 2C 2 (i + O + «(io + ttO) 



(20) 



The function -B(£) is determined by an expression The quantity a xxx has a static (u> — ► 0) limit: 

<WC, C) = -CS(C) 2 ^ - 8(4(1, C) + 2 (-1 + 6C 2 ) x 

m = I (-l+^ 2 )(-7r + 2^(l,0) $(2 C ) + 4 C( i _ 2 C 2 )$(3, C) + 2C 2 (C 2 - 1)4(4, C) J (22) 

^ 2 7r + C-7rC 2 + $(l,Oe(2$ 2 -3)' 1 J / 
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If additionally £ — > 0, 



If the frequency goes to infinity, 



<WC,C)« J-^Ciog(C/2)-C 



(4 + 3tt 2 
12tt 



... (23) 



The case of a clean sample gives the limit 1/6. The pos- 
itive slope of the function a xxx {(,, change to negative 
at very low numerical value of £ ~ 2exp(— 1 — 37r 2 /4). 
This behavior is plotted on the inset in Fig. 2. 
If C-»oo, 



+ ... 



(24) 



F(0 = 



4 ' 



1 

if 
if 



f(0, 



(25) 



( — > oo. 



The function -F(C) has the maximum equal to 0.0717 at 
C = 0.455. For finite u> and £ — > a« — ► 0. 
The static limit of is given by 



a xyy ((, C) = - 2tt 2 C(C 2 - I) 2 + tt (1 + 2C 4 - 6C 2 ) + 2 (C + 2( 5 - 2( 3 ) + C)(8C 6 - 16C 4 - 6C 3 + HC 2 ) 



4C 3 (2C 4 -4C 2 + 3) $ 2 (1,C) 



2 (C 2 - 1) (tt + C - < 2 + C*(l, C)(2C 2 - 3)) 



(26) 



which yields 



^y(C,C)«-2 + ^( 31 °g|+ 3 + 2 ^ 2 )+- if C^O, 

_^ ~ | (27) 

The function & xyy (C,, Q (similarly to 

&XXX 

(C, £)) has sin- 
gularity at ( = and changes the sign of slope at very 
low (. 

For large £ we have 



if C -> o, 



/ 1 4- 2fi 2 \ 

4 w ( C + ffl.0»-^( 5i ^- n |-l). <»> 



a xyy (C + i^,C) 



07T 

24C 5 



if £ — > oo. (28) The high-frequency behavior of <5 XOT/ (£ + iO, C) i s 
I 



Q: 



C 2 (l - C 2 ) 3tt - 4C - 2tt 2 C - 2ir( 2 + 2tt 2 C 3 + (8ir( 2 - tt + 4£ 3 - 8ttC 4 ) $(1, C) + 8£ 3 (C 2 - l) $(1, C) 2 



xyy 



20 2 



tt(£ 2 -!)-£ + £(3 -2£ 2 )$(1,£) 



(30) 



r 



with asymptotics 

, c 2 



and 



24£0 2 



In the static limit the component a yxy can be presented 
in the form 



«x W «-^2(3 + log|)... for C«l, (31) 



+ ... for £7>£»1. (32) 
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a yxy (Cj C) 



6C 2 (C 2 - i) 



2(C - 1)(2C 4 -1-C + 3C*(1, 0) " 3C ((C ~ 1)(16C - 11 



8C 4 + 4<(C 2 - l) 2 ) + (15C 2 - 6 - 20C 4 + 8C 6 )$(1, C))B(0 



(33) 



For small or large values of £ we have 
1 



OtyxyiCi 



3C 2 
n 



if C « 1, 
if C>1- 



For large f2 



ft 3 



(34) 



(35) 



Fi(C) = --( 20- 48C 2 + ^ 2 (1 - 4C 2 ) + 6ttC (-1 + 4( 2 ) 



a /a„ 

xxx 0.16670 



0.15- 



0.10- 



0.05- 



0.00 




0010 



1.5 £ 2.0 



+2 (26C 2 - 4 - 24C 4 + ttC (4( 2 - 3)) $(1, C)J , (36) 

M0 = ^ 2 (3C - 12C 3 ) - 96C (C 2 - 4) 
+2tt(25 - 88C 2 + 24C 4 ) 

-2C (144 + 9ttC - 224C 2 - 12ttC 3 + 48C 4 ) *(1, C))- (37) 

For small and large £ this gives 

■2 



n<-U'.„ lv i ^i(-5-^-21n(^)), if C^0, 



n 2y - 4 

7T 



FIG. 2: The dependence of a m on the parameter £ = t/tj 
for different frequencies f2 = ojt = 1,0.3,0.1,0.03,0 (marked 
on curves); Qo = e 3 /(4tt 2 vfN 2 D 2 ). Insert: the dependence 
of a xxx for small values of the parameter ( at = 0. 



if £ — ► 00. 



(38) 



1 ,25tt 



Im(5^)«^3(^+C(16 + y + 121n^)), if C - 0, 



(39) 



a /a 5 

xyy wo 




NUMERICAL RESULTS AND DISCUSSION 

The Figures f2!5l represent all components of &ijk cal- 
culated according Equations H18I2()|I at T = versus pa- 
rameter £ = r/rj for different frequencies. These depen- 
dencies can be treated as the dependencies of a%jk on the 
rate of impurity scattering. The sign of coefficient a xxx 
is positive, while the other coefficients change sign. In 
accord with found asymptotics all components tend to 
zero for £ — > 00 and exhibit non-analytical behavior at 

C = o. 

Asymptotic behavior a oc C 3 (or a rx t 1 ) at large £ 
follows from the odd dependence of the current on the 



FIG. 3: The dependence of &xyy on £ for the different param- 
eters f2 (marked on curves). 



asymmetric scattering on the cuts that results in the 
proportionality of the current to the scattering rate on 
cuts for low their concentration. This asymptotics cor- 
responds to case of weak asymmetric scattering usually 
considered in the theory of PGE. 

The value aijfe(0,0) depends on the order of limit 
Cij — > 0, 1/Tj — > 0: for example, if first to — > then 
1/tj — > a xxx — > 1/6, else a xxx — > 0. Such behavior re- 
sults from the absence of relaxation of electrons moving 
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FIG. 4: The dependence of linear photogalvanic coefficient 
Re(a yxy ) on £ for the same parameters as in the Figure 2. 



FIG. 6: Asymptotics of photogalvanic coefficient ot yxy = 
Ji(£)/f2 2 +iJ2(C)/^ 3 for large £7, according to the Equations 




FIG. 5: The dependence of circular photogalvanic coefficient 
Im^y^) for the same parameters il as in the Figure 2 



along y axis. The state with p x = plays role of a drain 
for electrons. These electrons do not participate in the 
transport along x axis, but due to absence of relaxation 
they accumulate in the state p x = 0; this suppresses the 
distribution function with finite p x , and j x — > 0. On the 
contrary, the transport along y axis diverges due to the 
same reason. 

For linear polarized electric field the signs of cur- 
rent components depend on the direction of polariza- 
tion. Physically, this can be explained by the effective 
increase of the mean-squared component of electron mo- 
mentum along field and subsequent increase (decrease) 
of scattering on the cuts. Let electron with a momentum 
p = (=tp, 0) impacts with a cut. The change of momen- 
tum are equal to — 2p for an electron with the momentum 
p = (p, 0) and (1 + 2/n)p for an electron with the oppo- 
site momentum, respectively. In equilibrium this change 
is compensated by the contributions of other electrons. 



But the increase of the mean-squared component of elec- 
tron momentum along field gives the finite positive con- 
tribution to the total current. The acceleration of an 
electron by the field in y direction increases y component 
of the momentum and produces the opposite direction of 

Jx- 

The values of coefficients a yxy are essentially larger 
than a xxx and a xyy . It is a consequence of the fact that 
the motion along y direction is collisionless unless the 
impurity scattering is taken into account. Obviously, the 
difference is more pronounced at low £ . 

It is useful to calculate the possible maximal PGE co- 
efficient. From said above it follows that this maximum 
is achieved for the component Re(a yxy ) at low frequen- 
cies a; -C 1/Tj and for clean material Ti r. In this case 
we have 

I / El \ ^ 
h = ~7^- n e ev F\ -J sin26>, (40) 

487T \ e F J 

where n e is the electron concentration, = VpTi, 9 is 
the angle of electric field with respect to the x direction. 
The estimation gives the value j y ~ 10~ 5 A/cm for U = 
10 _3 cto, rii = 10 12 cm~ 2 , vp ~ I0 7 cm/s. 

The special case is the limit of zero impurity concentra- 
tion. As we have found, in this case for finite frequency 
the current in x direction becomes equal to zero and in- 
finity for y direction. From said above it is evident that 
the inelastic collisions become essential: they limit this 
non-analytical behavior. The careful examination of this 
question is beyond the paper. 

The special case is the zero-frequency limit. As we 
have found, in this case for finite impurity scattering the 
current in x direction becomes equal to zero and infin- 
ity for y direction. From said above it is evident that 
the inelastic collisions become essential: they limit this 
non-analytical behavior. The careful examination of this 
question is beyond the paper. 
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We have restricted ourselves by the classical considera- 
tion only. Two quantum factors have not been taken into 
account: comparability of hw with characteristic electron 
energies, the Fermi energy and the temperature, and the 
quantum corrections to conductivity. It should be em- 
phasized that in used approximation the PGE in the 
metal case does not depend on the temperature in the 
low temperature limit. This is because the momentum 
relaxation in metal is temperature-independent (unlike 
energy relaxation). Despite the involvement of energy 
relaxation into the control on the distribution function, 
the current in the second order in electric field is de- 
termined by the contributions which do not depend on 
energy relaxation. 

The quantum corrections suppress the low tempera- 
ture transport and should lead to the tem pera ture de- 
pendence of the effect. Both modelling of |2fiL l27j and 
the present theory neglect quantum corrections. There- 
upon, the temperature dependence of steady current in 
[261 127| is not clear. The approach of these papers used 
a friction force essentially depending on the excess elec- 
tron energy. Possibly, it is the reason of the temperature 
dependence. It should be mentioned that the quantum 
corrections are of less importance in high-mobility sam- 
ples utilized usually in experiments with antidots. 

COMPARISON WITH NUMERICAL 
SIMULATIONS 

In this section we compare our results with the sim- 
ulations of the flow velocity of a particle colliding with 
semidisks. We consider the model used here as a sim- 
plification of the semidisks model of Ref. 1281. The es- 
sential differences between approaches of Ref. |28( and the 
present paper are: periodic/random distribution of asym- 
metric scatterers, deterministic/chaotic character of elec- 
tron scattering on antidots and motion between them, 
and absence/presence of the gas approximation. How- 
ever, the case of low density of semidisks and weak fric- 
tion corresponds to the gas approximation because the 
scattering randomizes the motion well enough, and the 
results of both approaches should be in accord with each 
other. 

The Fig. 6 in |28j] shows approximate quadratic depen- 
dence of the flow velocity on the alternating field what 
agrees with the second-order in field approximation used 
in the present paper. 

The dependence Fig. 7 from demonstrates the drop 
of the flow velocity with the growth of the semidisk 
size. This behavior qualitatively corresponds to the drop 
ctxxx — > for t — ► (Fig. 2 of present paper). At the 
same time there is no drop in this figure for large r — * oo 
following from the present paper. 

The dependence of the flow velocity on the distance 
between semidisk centers (Fig. 8 from 28]) can be treated 



as rescaled dependence of a xxx on the mean free time r 
which in the present case vanishes both for r — > and 
t — ► oo (Fig. 2 of the present paper). 

The dependence of the flow velocity (Fig. 9 in |2^) on 
the impurity scattering also has a drop for Tj — ► 0, like 
expected from the Ea. l|24|l . but there is no drop in this 
figure for large r%. 

Hence, there is a partial qualitative accordance be- 
tween the present results and computer simulations of 
pjj . The origin of discrepancy needs additional study. 



CONCLUSIONS 

The considered system of oriented asymmetric 
scatterers-cuts has a xxx , a xyy , and a yxy non-zero compo- 
nents of photogalvanic tensor. The linear photogalvanic 
effect is determined by a xxx , ct xyy , and Ke(a yxy ). The 
circular-polarized illumination, causes the response of the 
y component of current only, determined by Im(a yxy ). 
The static limit of the current in impurity-free system 
is ambiguous, depending on value of the product of fre- 
quency to the impurity mean-free time. The x compo- 
nent of the current is limited and in the impurity-free 
system tends to zero, while the y-component tends to in- 
finity. This is explained by the accumulation of electrons 
in the state with zero rr-component of momentum. 
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